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ABSTRACT 


--  It  is  found  that  it  is  possible  to  propagate  a  horizontally  polarized 
(SH)  wave  without  dispersion  through  an  elastic,  periodically-layered  com¬ 
posite.  The  nondispersive  property  of  the  wave  is  due  to  the  fact  that  at 
each  interface  the  angle  of  incidence  is  such  that  the  wave  is  totally  trans¬ 
mitted  without  reflection.  In  optics  such  an  angle  is  referred  to  as  the 
Brewster  angle.  It  is  determined  that  this  particular  case  is  contained  as 
a  special  solution  of  the  general  dispersion  equation  for  SH  waves,  which  has 
not  been  noticed  before. 
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1.  INTRODUCTION 

For  horizontally  polarized  shear  waves  impinging  upon  an  interface  be¬ 
tween  two  media  with  different  properties,  it  is  well  known  that  in  general 
there  will  be  both  a  transmitted  and  a  reflected  wave.  For  a  special  angle  of 
incidence  the  amplitude  of  the  reflected  wave  is  known  to  be  zero  and  only 
the  transmitted  wave  will  remain.  In  optics  this  angle  is  referred  to  as  the 
Brewster  angle.  Under  these  very  special  conditions  it  can  be  expected  that 
this  same  reflection-free  phenomenon  must  exist  in  a  bi-laminated  composite. 
This  is  a  special  case  which  has  been  overlooked  in  the  general  treatment  of 
wave  propagation  in  anti-plane  strain  [1].  It  is  the  aim  of  this  paper  to  ex¬ 
amine  in  detail  this  particular  phenomenon  of  nom  spersive  wave  propagation 
in  a  layered  medium. 
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BREWSTER  ANGLE 


We  begin  by  examining  the  case  of  an  infinite  SH  wave  train  impinging  at 
an  angle  0  upon  the  interface  of  two  semi-infinite  elastic  homogeneous  media 
of  different  properties.  In  general,  a  portion  of  the  wave  will  be  reflected 
and  a  portion  will  be  transmitted  (see  Fig.  la).  For  a  particular  angle  of  in¬ 
cidence  the  wave  will  be  totally  transmitted  (Fig.  lb).  The  angle  8g  is 
called  the  Brewster  angle  and  is  given  by  an  expression  of  the  form  [2] 


sin  eB  = 


(1) 


where  y,p  and  p',p'  are  the  shear  moduli  and  density  in  the  two  media, 
respectively.  The  feature  of  no  reflection  can  be  generalized  from  the  case 
of  two  semi-infinite  media  in  contact  to  the  case  of  an  infinite  laminated 
medium  (Fig.  2).  Here  the  condition  of  no  reflection  occurs  at  each  inter¬ 
face.  As  a  result  we  can  state  that  the  laminated  composite  is  capable  of 
supporting  such  reflection-free  waves.  An  immediate  observation  is  that  these 
waves  are  independent  of  the  frequency,  i.e.  "nondispersive" .  In  frequency- 
wavenumber  space,  waves  propagating  in  such  a  fashion  are  represented  by  a 
straight  line.  A  discussion  of  this  phenomenon  for  antiplane  strain  wave 
propagation  in  an  infinite  periodically  layered  composite  has  not  been  noted 
before.  In  the  next  section  we  will  take  a  closer  look  at  the  dispersion 
relation  for  the  infinite  laminated  composite  and  examine  this  noteworthy  dis¬ 
persion-free  propagation  phenomenon. 
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3.  DISPERSION  IN  AN  INFINITE  LAMINATED  MEDIUM 

We  consider  horizontally  polarized  harmonic  shear  waves  (SH  waves)  prop¬ 
agating  through  a  periodically  layered  elastic  body  of  unbounded  extent.  The 
union  of  any  two  contiguous  layers  in  the  body  constitutes  a  unit  oell,  and 
this  unit  cell  is  invariant  under  a  lattice  translation  along  the  positive 
and  negative  y-axes  (Fig.  2).  Each  of  the  two  layers  in  the  unit  cell  are 
assumed  to  be  homogeneous,  isotropic,  and  perfectly  bonded  to  the  adjoining 
layers.  The  two  lamellae  of  a  typical  unit  cell  have  elastic  constants 
[(A,p)  ;  (X  * » u  * )  ]  «  thicknesses  (2h;2h')  ,  and  densities  (p  ;  p  * )  ,  respectively. 

Let  u  ,  v  and  w  be  the  three  Cartesian  components  of  the  displace¬ 
ment  vector  in  the  x,  y  and  z-directions ,  respectively.  For  antiplane  motion, 
we  take 

v  =  w  =  0  and  u  =  u(y,z;t)  •  (2) 

For  the  layers  with  unprimed  constants  in  the  Nth  unit  cell,  the  equation  of 
motion  is 


V 


("h  -  yN  -  h)  *  c  >  0 


(3) 


where  yN  is  a  local  coordinate  with  its  origin  at  the  midplane  of  the  layer 
(Fig.  2).  For  plane  waves  traveling  in  the  positive  z-direction,  we  assume 
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where  k  is  the  wave  number  in  the  z-direction  and  w  is  the  circular  fre- 
z 

quency  in  radians  per  unit  of  time.  Substitution  into  (3)  yields  the  ordinary 
differential  equation 


+  (‘r-  -  kz V  ' 0  •  s  i 
<‘>'8  '  / 


whose  solution  may  be  written  in  the  form 


f(VN)  =  Cx  exp  /  ■ -—U-  C2  exp 


(-?) 


o  =  Vft2  -  c2  , 


and  ft  and  C  are,  respectively,  the  nondimens ional 


frequency  and  wave  number  defined  by 


2hw  2h 

a  .  — —  ;  c  -  -  \ 


The  displacement  component  now  takes  the  form 


<  h)  ,  t  >  0 


(8) 
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In  a  similar  manner,  the  displacement  component  of  the  layer  with  primed  con¬ 
stants  in  the  Nth  unit  cell  can  be  written  as 


x  exp 


/  i’o  y»  \  /  \  1 
IK(‘Z  -^f“C)]  •  ‘-h'iyNih'> 


t  >  o 


(9) 


where  a'  =  V and  =  (pp'/y'p)  . 

The  layers  are  considered  to  be  perfectly  bonded,  so  continuity  of  dis¬ 
placement  and  traction  must  be  enforced  at  the  layer  interfaces.  The  conti¬ 
nuity  conditions  between  the  two  layers  in  the  Nth  unit  cell  and  the  adjoining 
layer  in  the  (N  +  l)th  unit  cell  are 


u(-h,zN;t)  =  u'(h',z';t)  , 


°xy(-h*ZN;t)  =  0xy(h'’zN;t)  ’ 


u(h,zN;t)  =  u'C-h'.z'^jt)  , 


"  VZN+1  <  “  >  t  >  0 


axy(h,zN;t)  =  °xy(-h'»zj5+1;t)  • 


(10) 


It  may  now  be  noted  that  the  equation  of  motion  for  the  laminated  body 
B  takes  the  form  of  a  partial  differential  equation  with  coefficients  which 
are  periodic  in  the  y-direction  with  period  d  ■*  2(h  +  h')  ,  This  follows 
from  the  fact  that  the  coefficients  of  the  equation  of  motion  (the  elastic 
moduli  and  mass  density)  are  piecewise  constant  in  each  layer  and  have  periodic 
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variation  from  cell  to  cell.  It  is  well  known  from  the  one-dimensional  theory 
of  Floquet  [4]  or  the  three-dimensional  theory  of  Bloch  that  differential 
equations  with  periodic  coefficients  admit  solutions  of  the  form 


u(y,z;t)  =  g(y,z;t)  elkyy 


(11) 


where  g(y,z;t)  has  the  same  periodicity  as  the  coefficients  of  the  differ¬ 
ential  equation,  i.e., 


g(y  +  d,z;t)  =  g(y,z;t)  . 


(12) 


From  (11)  it  is  seen  that  has  the  character  of  a  wave  number;  it  is  said 

to  be  the  wave  number  of  the  Floquet  wave.  Equations  (11)  and  (12)  lead 
immediately  to  the  quasi-periodic  recurrence  relation 


u(y  +  d,z;t)  =  u(y,z;t)  elkyd 


(13) 


The  Floquet  solution  (11)  has  several  other  interesting  properties.  In 
particular,  it  is  not  difficult  to  show  that  the  wave  number  k^  is  uniquely 
determined  only  to  within  an  integer  multiple  of  2ir/d  .  This  fact,  its 
physical  significance,  and  other  properties  of  the  Floquet  solution  are  dis¬ 
cussed  in  some  detail  by  Lee  [3]. 

Using  (13),  the  last  two  equations  of  (10)  may  be  written  as 


u(h,zN;t)  =  u  (-h',z^;t)  eikyd 

a  (h,z  ;t)  =  o'  (-h',z';t)  eikyd 
xy  N  xy  N 


-®>  <  z. 


t  >  0 


(14) 
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which  together  with  the  first  two  equations  of  (10)  now  represent  a  set  of 
displacement  and  traction  continuity  conditions  written  in  terras  of  the  dis¬ 
placements  and  stresses  at  the  interfaces  of  the  Nth  unit  cell.  If  these  are 
satisfied,  the  quasi-periodicity  condition  (13)  will  insure  continuity  of  dis¬ 
placement  and  traction  across  the  interfaces  of  every  unit  cell  of  the  lami¬ 
nated  medium. 

Substituting  (8)  and  (9) ,  and  the  isotropic  stress-strain  relation 

o  =  u(3u/3y)  into  (14)  and  the  first  two  equations  of  (10),  we  obtain  a 
xy 

set  of  four  homogeneous  equations.  For  a  nontrivial  solution,  the  determinant 
of  the  matrix  of  coefficients  must  vanish,  yielding  the  dispersion  equation 


e* 

-ee*' 

-e‘e*' 

-<s> 

yae 

<P 

-yae 

..  E  Ip' 

-a  e 

e* 

o’* 

_ed-f.4>') 

(ip+e<p  *) 

-e 

=  0 

(15) 

<P 

yae 

-4, 

-yae 

-a'e(,H:',r) 

>  (4+e<t’ ') 

a  e 

ra/2  , 

<t>'  =  iita'/2  , 

=  iiT(l+e)n  » 

e  =  h'/h  ,  y  = 

v/u'  , 

and 

This 

equation  may 

be  written  in 

the  convenient  form 

=  4yaa 

cos  irrid  +  e 

)  +  (ya  -  a') ^ 

cos  ir(a  -  ea') 

-  (ya  +  a')* 

cos  Tr(a  +  ea') 

=  0  . 

(16) 

Since  this  equation  remains  unaffected  by  a  change  in  sign  in  q  or  5  ,  both 
these  quantities  may  be  assumed  positive.  However,  if  either  one  or  both  of 
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a  and  a'  is  imaginary,  (16)  must  be  appropriately  interpreted  in  terms  of 
hyperbolic  functions. 

It  may  be  noted  that,  since  (16)  relates  a  nondimensional  frequency  to 
two  nondimensional  wave  numbers,  the  roots  of  this  equation  define  a  surface 
in  frequency-wave  number  space,  which  is  called  the  dispersion  surface.  A 
qualitative  sketch  of  this  surface,  drawn  on  an  extended  zone  scheme  [3],  is 
shown  in  Fig.  3.  It  can  be  seen  that  the  surface  is, in  general, discontinuous 
at  n  =  n/(l  +  e)  ,  n  =  1,2,....  These  planes  of  discontinuity  divide  the  sur¬ 
face  into  Brillouin  zones,  the  first  three  of  which  are  shown  in  the  sketch. 

The  geometric  and  material  parameters  used  here  and  throughout  the  remainder 

2 

of  this  paper  are  y  =  0.02  ,  a  =  0.06  ,  and  c  =  4  . 

In  the  past  the  physical  interpretation  of  the  wave  numbers  n  and  c  has 
been  quite  deficient.  Until  now  no  simple  example  was  known  which  would  give 
insight  into  the  nature  of  these  wave  numbers. 

A  description  of  the  nondispersive  propagation  in  terms  of  the  wave- 
numbers  n  and  c  illustrates  a  point  which  is,in  general , not  sufficiently 
appreciated.  In  some  papers  (e.g.  [5]),  the  wavenumbers  n  and  £  were 

treated  as  vector  components  of  a  single  wave  number  measured  normal  to  a  plane 
wave  front.  It  is  clear  from  Fig.  5  that,  in  general,  the  wave  forms  of  the 
type  appropriate  to  layered  composites  are  not  plane  waves,  except  in  the  case 
C  “  0  .  Intuitively,  some  average  propagation  direction  should  be  described 
by  the  relation  C/n  .  Our  intuition  fails  us  if  we  guess  that  this  direction 
may  be  defined  by  connecting  any  two  corresponding  points  in  adjacent  unit 
cells  in  the  wave  vector  representation  of  Fig.  4a.  Based  on  the  previous 
numerical  values  chosen,  we  can  calculate  an  effective  propagation  direction 
as  defined  by  the  line  connecting  two  corresponding  points  (A,B)  in  adjacent 


unit  cells  from  the  wave  vector  representation  (see  Fig.  4b) .  For  the  para¬ 
meters  chosen, an  angle  of~6°  is  obtained  for  0  . 

By  careful  inspection  of  the  dispersion  surface  we  find  that  the  only 
possible  direction  for  a  straight  line  in  frequency-wave  number  space  which 
lies  entirely  on  the  surface  is  illustrated  by  the  straight  line  OABC  in 
Fig.  2.  It  can  be  shown  that  the  equation  of  this  line  satisfies  the  disper¬ 
sion  equation.  Again,  for  the  particular  parameters  chosen  we  find  for  a 
propagation  direction  given  by  4/ n  a  value  of  ~41°.  Thus  our  initial  guess 
for  the  propagation  direction  based  on  the  wave  vector  representation  is  not 
correct . 

A  dual  representation  to  the  wave  vector  is  by  means  of  the  wave  front. 

In  an  infinite  homogeneous  isotropic  elastic  space  the  wave  vector  is  normal 
to  the  wave  front.  For  the  laminated  composite  this  is  true  within  each  layer, 
but  not  in  an  average  sense.  Fig.  5a  shows  the  wave  front  corresponding  to 
the  wave  vectors  of  Fig.  4a.  Fig.  5b  shows  the  "averaged"  wave  front.  The 
calculated  angle  0  for  the  assumed  properties  was  found  to  be  identical  to 
the  value  obtained  from  the  dispersion  surface.  This  can  be  shown  mathemat¬ 
ically  by  considering  two  such  wave  fronts  as  shown  in  the  next  section. 

This  shows  quite  clearly  the  manner  in  which  n  and  4  can  be  regarded 
as  vector  components  of  some  "equivalent"  wave  front.  This  illustrative  ex¬ 
ample  derived  from  this  nondispersive  propagation  direction  appears  to  be  the 
only  case  where  wave  fronts  can  be  drawn  easily. 

An  unexpected  sidelight  of  this  analysis  deals  with  the  concept  of  what 
is  meant  by  a  nondispersive  wave.  We  previously  stated  that  a  straight  line 
in  frequency-wave  number  space  represents  nondispersive  wave  propagation.  This 


could  be  interpreted  as  meaning  the  group  velocity,  as  defined  by 
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3H 

3fi 


remains  constant.  This  is  in  addition  to  the  phase  velocity 


where  e^  and  e^  are  unit  vectors  perpendicular  and  parallel  to  the  layers, 
respectively,  remaining  constant.  For  the  present  case  the  line  OABC  in 
Fig.  3  certainly  satisfies  the  second  criterion.  But  the  group  velocity  is 
only  constant  if  we  restrict  our  attention  to  the  line  OABC  .  In  an  experi¬ 
ment  in  which  the  .  utio  of  5/n  is  fixed  and  we  move  only  along  the  line 
OABC  ,  we  would  indeed  find  that  the  system  is  nondispersive .  If,  instead, 
one  is  permitted  to  sample  in  other  propagation  directions,  we  would  find 
that  the  group  velocity  is  not  constant  along  the  line  OABC  and  the  disper¬ 


sive-free  nature  is  lost. 


1 


4.  SPECIAL  SOLUTION  OF  THE  GENERAL  DISPERSION  EQUATION 

In  this  section  we  will  show  that  there  exists  a  straight  line  in  fre¬ 
quency-wave  number  space  which  satisfies  the  dispersion  equation  identically. 
Also  it  will  be  shown  that  the  ratio  of  the  wave  numbers  for  the  effective 
plane  wave  is  the  same  as  the  ratio  of  £  to  n  found  previously. 

We  first  show  that  the  equation  of  a  line  through  two  points  ,  £  , 

tig)  and  ,  r^)  in  frequency-wave  number  space  can  be  written  in  the 

following  form: 


For  convenience  we  choose  ,  nQ)  =(0,0,0)  and  let  (£2^  ,  £  , 
n1>  be  the  coordinates  of  the  point  A  in  Fig.  3.  In  this  case  tht  equation 
of  the  line  through  these  two  points  can  be  written  as: 


(18) 


We  can  rewrite  the  dispersion  equation  in  the  following  form: 
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From  (18)  we  find  for  0/£ 


(20) 


Substituting  into  (19)  and  dividing  by 


2 


gives,  after  rearranging. 


But  from  (18)  we  see  that 


n(l  +  e) 


(1  +  ye)? 


(22) 


In  this  case  the  dispersion  equation  is  satisfied  identically  for  the  line 
defined  by  (18). 
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The  only  remaining  question  concerns  the  relationship  between  the  wave 
number  components  of  the  effective  plane  wave  and  the  wave  numbers  5  and  n  • 
To  clarify  this  point  we  consider  two  consecutive  wave  fronts  defined 
by  BIF  and  EDH  (see  Fig.  6).  The  straight  lines  BF  and  EH  are  the  effective 
plane  wave  fronts.  Using  some  definitions  from  wave  mechanics  and  trigonom¬ 
etry,  we  will  now  derive  the  relation  between  k  and  k  (the  effective  wave 

x  y 

numbers  in  the  horizontal  and  vertical  direction,  respectively)  and  the  wave 
number  £  and  n  .  We  note  first  that 

\  _  EF  =  _ 2(h  +  hQ _ 

ky  FH  2 (h  cot  B  +  h'  cot  B")  * 

This  can  be  rewritten  as 


k 

x 

k 

y 


cot  e ' 


1  +  e 


+ 


(2*) 


i<  ^  ^  ^  \< 

where  z  =  (pp)  ,  z'  =  (p 'v')  are  the  impedances,  and 

-  u 

Cp  =  (v'/p')  are  the  phase  velocities.  The  expression 
be  rewritten  as 


cp  =  (w/p)1* 
in  parentheses  can 


(25) 


Substituting  (25)  into  (24)  gives 
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1  +  e  _  y(l  +  e)  tan  B' 


ky  cot  +  t) 


( 1  +  ye) 

Using  Snell's  Law  we  can  obtain,  after  some  manipulation 


(26) 


-  1 


tan  8'  = 


c 

p 


(27) 


This  can  be  rewritten  as 


-  1’ 
2  1 


tan  &' 


1  -  o' 


I 
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(28) 


Thus  we  obtain  after  substituting  into  (26)  and  rearranging 


(1  +  e) 


1  +  ye 


(29) 


If  we  now  return  to  equation  (18),  we  note  that  the  expression  for  4/n  is 
identical  to  equation  (29).  Thus  finally  we  have 


(30) 
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FIGURE  CAPTIONS 


Fig.  la  SH  wave  impinging  upon  interface 

Fig.  lb  SH  wave  impinging  at  Brewster  angle 

Fig.  2  Geometry  of  layered  composite 

Fig.  3  Dispersion  surface  for  laminated  medium 

Fig.  4a  Wave  vector  representation 

Fig.  4b  Effective  wave  vector 

Fig.  5a  Wave  front  representation 

Fig.  5b  Effective  wave  front 

Fig.  6  Detail  of  wave  front 
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